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ABSTRACT 

A general model is introduced to describe a wave-envelope system for the situation when the linear 
dispersion relation has three branches, which in the absence of any coupling terms between these branches, 
would intersect pair-wise in three nearly-coincident points. The system contains two waves with a strong 
linear coupling between them, to which a third wave is then coupled. This model has two gaps in its 
linear spectrum. As is typical for wave-envelope systems, the model also contains a set of cubic nonlinear 
terms. Realizations of this model can be made in terms of temporal or spatial evolution of optical fields 
in, respectively, either a planar waveguide, or a bulk-layered medium resembling a photonic-crystal fiber, 
which carry a triple spatial Bragg grating. Another physical system described by the same general model 
is a set of three internal wave modes in a density-stratified fluid, whose phase speeds come into close 
coincidence for a certain wavenumber. A nonlinear analysis is performed for zero-velocity solitons, that 
is, they have zero velocity in the reference frame in which the third wave has zero group velocity. If 
one may disregard the self-phase modulation (SPM) term in the equation for the third wave, we find an 
analytical solution which shows that there simultaneously exist two different families of solitons: regular 
ones, which may be regarded as a smooth deformation of the usual gap solitons in a two-wave system, 
and cuspons, which have flnite amplitude and energy, but a singularity in the first derivative at their 
center. Even in the limit when the linear coupling of the third wave to the first two nearly vanishes, 
the soliton family remains drastically different from that in the uncoupled system; in this limit, regular 
solitons whose amplitude exceeds a certain critical value are replaced by peakons. While the regular 
solitons, cuspons, and peakons are found in an exact analytical form, their stability is tested numerically, 
which shows that they all may be stable. If the SPM terms are retained, we find that there may again 
simultaneously exist two different families of generic stable soliton solutions, namely, regular ones and 
peakons. Direct simulations show that both types of solitons are stable in this case. 

PACS numbers: 05.45.Yv; 42.65.Tg; 42.81.Dp; 47.55.Hd 
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1 Introduction 



1.1 The model system 

Gap solitons (GS) is a common name for solitary waves in nonlinear models which feature one or more 
Q gaps in their linear spectrum, see review papers Q and respectively, for these two cases. A soliton 
may exist if its frequency belongs to the gap, as then it does not decay into linear waves. 

Gaps in the linear spectrum are a generic phenomenon in two- or multi-component systems, as in- 
tersection of dispersion curves belonging to different components is generically prevented by a linear 
coupling between the components. Excluding cases when the system's linear spectrum is unstable (which 
are possible in a fluid dynamics application ||^), the intersection avoidance alters the spectrum so that a 
gap opens in place of the intersection. Approximating the two dispersion curves, that would intersect in 
the absence of the coupling, by straight lines, and assuming a generic cubic [x^^"*] nonlinearity, one arrives 
at a generalized massive Thirring model (GMTM), which has a family of exact GS solutions that com- 
pletely fill the gap 1^. The model has a direct application to nonlinear optics, describing co- propagation 
of forward- and backward-traveling electromagnetic waves in a fiber with a resonant Bragg grating (BG) . 
Gap solitons, first predicted theoretically, were observed in experiments with light pulses launched into a 
short piece of the BG-equipped fiber Q (in fact, optical solitons that were first observed in the BG fiber 
were, strictly speaking, not of the GS type, but more general ones, whose central frequency did not 
belong to the fiber's bandgap). 

GS are known not only in optics but also in other physical settings, for instance, in density-stratified 
fluid flows, where dispersion curves pertaining to two different internal-wave modes often exhibit near- 
intersections. Again taking into regard weak nonlinearity, one can predict the occurrence of GS in 
density-stratified fluids . 

In this work, we aim to consider GS that may exist in a generic situation of the next type, when the 
underlying system contains three wave components, and the corresponding dispersion curves intersect 
at three nearly-coincident points, unless linear coupling terms are taken into account. Situations of this 
type can readily occur in the above-mentioned density-stratified fiuid fiows, since tuning of two suitable 
external parameters can often lead to a near-coincidence in the linear phase speeds of three independent 
internal wave modes, for certain wavenumbers (see |^). Indeed, similar considerations can be applied 
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to many other physical systems. As for the usual GS systems, the generic wave-envelope model can be 
expected to contain cubic nonlinear terms. 

In optics, a ^'^'^^ -nonlinear model with three linearly coupled waves is possible too, in terms of either 
temporary evolution of fields in a planar nonlinear waveguide equipped with a triple BG in the form of 
three systems of parallel scores, or spatial evolution of stationary fields in a bulk waveguide with a similar 
triple BG consisting of three systems of parallel interfaces between layers. The latter realization seems 
natural enough, as it strongly resembles photonic-crystal fibers, which have recently attracted a great 
deal of interest . Note that the former version of the model is a generalization of a three- wave model 
for a x'^'-nonlinear planar waveguide with an ordinary BG, which was introduced in [p^ . Both versions 
of the proposed model are illustrated in Fig. 1, where the periodic lattice shows the triple BG. In the 
case of the temporal evolution. Fig. 1 displays the planar waveguide, while in the case of the spatial 
evolution, it is a transverse cross-section of the bulk waveguide. 

We stress that the lattice in Fig. 1 is not completely symmetric; although the triangular cells of the 
lattice are equilateral ones, the two diagonal sub-gratings are assumed, in the general case, to have the 
strength (contrast of the refractive index) smaller than the horizontal one. The bold triangles inscribed 
into the two triangular cells illustrate the resonant Bragg reflections that give rise to linear couplings 
between the waves. Then, neglecting intrinsic dispersion or diffraction of the waves in comparison with the 
strong artificial dispersion/diffraction induced by the Bragg reflections, normalized equations governing 
the spatial evolution of the three fields whose Poynting vectors are shown by three bold arrows in Fig. 
1, are 

/ dui dui 1 dui\ /, ,9 „, ,9 „, ,,s „ 

du2 , du2 1 9U2\ , , * , /I |2 , oi ,2 , ol \2\ 



i i-^ + ^-^) + + '^"S + (|W3P + 2|miP + 2|U2P) U3 = UJ0U3 . (3) 

Here, the evolution variable t is the proper time in the case of the temporal evolution in the planar 
waveguide, or the coordinate z in the direction perpendicular to the plane of the figure in the case of 
the spatial evolution in the bulk waveguide. In the latter case, the beam enters the medium through the 
plane z — and evolves along the coordinate z, that is why it plays the role of the evolution variable in 
the equations (the initial conditions necessary to launch a soliton will be discussed in more detail below). 
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The relative coefficients in front of the x- and y- derivative terms correspond to the geometry in Fig. 1, 
the coefficient of the walk-off in the x-direction in the first two equations being normalized to be ±1. The 
coefficient of the BG-induced linear conversion between the waves ui and U2 is normalized to be 1, while 
the parameter k (which is complex, in the most general case, but see the discussion below) accounts for 
the linear conversion between these waves and the third wave U3, and the usual ratio 1 : 2 between the 
coefficients of the self-phase modulation (SPM) and cross-phase modulation (XPM) is adopted. Lastly, 
(jjQ is a frequency/wavenumber mismatch between the third and the first two waves, which is caused by 
the above-mentioned asymmetry between the diagonal and horizontal subgratings, as well as by other 
reasons. 

As we mentioned above, the Bragg constant k in Eqs. (|l|) - (||), which couples the field M3 to the pair 
Ml. 2, is complex in the general case (note that the constant of the Bragg coupling between the fields ui 
and U2 might also be complex in its primary form, and making it equal to 1 in Eqs. (|^) and (^ involves 
opposite constant phase shifts of the fields ui and U2, which is why k and k* appear exactly as in Eqs. (|l|) 
- (|3|)). However, assuming that each score, the families of which constitute the triple grating (in the case 
of the temporal-domain evolution in the planar waveguide), gives rise to simple refiection described by 
the classical Fresnel formulas, it is easy to conclude that all the coupling constants are real and positive, 
provided that either the light is polarized orthogonally to the waveguide's plane, and the reflection takes 
place from a less optically dense material (i.e., the "score" is, literally, a shallow trough on the surface of 
the planar waveguide), or the light is polarized parallel to the plane of the waveguide, and the reflection 
is from a more optically dense material. Similarly, in the case when the same equations describe the 
spatial evolution of the optical flelds in the layered bulk medium, one may assume that either the light 
is polarized in the z -direction, and seams between the layers are filled with a material (for instance, air) 
which is optically less dense than the bulk medium, or the polarization is orthogonal to the z-axis (i.e., 
it is parallel to the plane of Fig. 1), and the material filling the inter- layer seams is optically denser than 
the host medium. 

In the present paper, we focus on this case, which was described above in detail for the realizations 
of the model in terms of both planar and bulk optical waveguides, and which corresponds to k real and 
positive in Eqs. (|^) - (p. Note, incidentally, that the case when k is real and negative can be reduced to 
the same case simply by reversing the sign in the definition of U3. 
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The model displayed in Fig. 1 may be further generalized by introducing an additional asymmetry, 
which will remove the equality between the horizontal side of the lattice's triangular cell and its diagonal 
sides. Then, the simultaneous fulfillment of the Bragg-reflection conditions for the waves 1*1^2 and U3 can 
be secured by making the waveguide anisotropic. However, such a generalization goes beyond the scope 
of this paper. 

For the physical realization of the model, Eqs. (|^) - (^) must be supplemented by initial conditions at 
t = in the case of the temporal evolution in the planar waveguide, or boundary conditions at z = in 
the case of the spatial evolution in the bulk medium. It is sufficient to assume that, at t = 0, a single wave 
component (for instance, U3) is launched into the waveguide. The linear-coupling terms in the equations 
will then start to generate the other components, and, if solitons that might exist in this model are stable 
(see below), they may self-trap from such an initial beam. 

Bearing in mind also the above-mentioned application to internal waves in stratified fluids, as well as 
similar realizations in other physical media, Eqs. (|^) - (^ may be naturally extended by introducing more 
general SPM and XPM coefhcients, as in applications other than nonlinear optics, the ratios between 
the XPM and SPM coefficients may be different from those adopted above. Thus the generalized system 
of equations takes the following form, in which we confine consideration to y-independent solutions (the 
consideration of possible three-dimensional solitons in the case of y-dependent fields is not an objective 
of this work) , 

- +U2 + KU3 + a (acTiluip + a|u2p + luaP) = , (4) 

j(^^ + + ui + KU3 + a (acri|u2p + + lusP) -"2 = 0, (5) 

du 

i-Q^ + K (ui + U2) + (cTslusI^ + a\uif + a\u2f) U3 = UJ0U3 . (6) 

where, in accord with the discussion above, we set k to be real and positive. 

The coefficients cri.s and a in Eqs. (^) - (^ are the generalized SPM and XPM coefficients, respectively. 
In particular, a is defined as a relative XPM coefficient between the first two waves and the third wave. In 
fact, the coefficients cti and both may be normalized to be ±1, unless they are equal to zero; however, 
it will be convenient to keep them as free parameters, see below (note that the SPM coefficients are always 
positive in the optical models, but in those describing density-stratified fluids they may have either sign). 
In optical models, all the coefficients a and cti^s are positive. However, in the models describing the 



internal waves in stratified fluids, there is no inherent restriction on their signs, and some of them may 
indeed be negative. 

The symmetry between the walk-off terms in Eqs. and is not really essential, and we will 
comment later on the more general case when these terms are generalized as follows: 

dui dui ^du2 ^ dui 

dx dx ' dx dx ' 

where ci and C2 are different, but have the same sign. As for Eq. (ref3), it is obvious that the walk-off 
term in this equation, if any, can always be eliminated by means of a straightforward transformation. 

We have kept only the most natural nonlinear SPM and XPM terms in Eqs. - (||), i.e., the terms 
of the same types as in the standard GMT model. Additional terms, including nonlinear corrections to 
the linear couplings (e.g., a term ~ |uipU2 in Eq. (^) may appear in more general models, such as a 
model of a deep (strong) BG ]Tot . 

Equations (|) - (|) conserve the norm, which has the physical meaning of energy in optics, 

/ + 00 
\Un{x)f dx, (8) 

the Hamiltonian, 

= -^grad ^ -^coupl ^" -^focusi (9) 



/• + CO 

■ffcoupl = [ulu2 + Kul{ui+U2)]dx + C.C., (11) 



— OO 



^4:1,2) 



and the momentum, which will not be used here. In these expressions, the asterisk and c.c. both stand 
for complex conjugation, -ffgiad, -ffcoupi, and iJfocus being the gradient, linear-coupling, and self-focusing 
parts of the Hamiltonian. To obtain the Eqs. (^ - (^from the Hamiltonian, the conjugate pairs of the 
variables are defined, in a standard fashion, as , ii„ . 

Our objective is to find various types of solitons existing in the generic three-wave system (^) - (^ 
and investigate their stability. Focusing first on the case (suggested by the analogy with GMTM) when 
the SPM term in Eq. (^ may be neglected (i.e., (T3 = 0), in section 3 we find a general family of 
zero- velocity solitons in an exact analytical form. We will demonstrate that they are of two drastically 



different types: regular GS, and cuspons, i.e., solitons with a cusp singularity at the center, in which 
the soliton amplitude is finite, but the derivative is infinite; further, the energy of the cuspons is finite. 
Cuspons are known to exist in degenerate models without linear terms (except for the evolution term such 
as du/dt), i.e., without a linear spectrum, a well-known example being the exactly integrable Camassa- 



Holm (CH) equation |17, (see also IQ). Our model resembles the CH one in the sense that both 
give rise to coexisting solutions in the form of regular solitons and cuspons. The cause for the existence 
of these singular solitons in our model is the fact that, looking for a zero-velocity soliton solution, one 
may eliminate the field U3 by means of an algebraic relation following, in this case, from Eq. (^). 
The subsequent substitution of that result into the first two equations (^ and (|^) produces a rational 
nonlinearity in them, the corresponding rational functions featuring a singularity at some (critical) value 
of the soliton's amplitude. If the amplitude of a regular-soliton solution is going to exceed the critical 
value, it actually cannot exist, and in the case when (73 = it is replaced by a cuspon, whose amplitude 
is exactly equal to the critical value. 

In the limit k — > 0, which corresponds to the vanishing linear coupling between the first two and third 
waves, the cuspon resembles a peakon, which is a finite-amplitude solitary wave with a jump of its first 
derivative at the center. Note that peakon solutions, coexisting with regular solitons (this property is 
shared by our model), are known in a slightly different (also integrable) version of the CH equation, see, 
e.g., Refs. 17, 2^, |l). We also note that soliton solutions with a discontinuity in the first derivative have 
been found in the BG model (which does contain a linear part) in the case where the grating parameter 
changes abruptly 

Then, we show that, when the SPM term is restored in Eq. (^) (i.e., 0-3 ^ 0; the presence or absence 
of the SPM terms ~ cri in Eqs. (jj) and (||) is not crucially important), the system supports a different 
set of soliton solutions. These are regular GS and, depending on the sign of certain parameters, a family 
of peakons, which, this time, appear as generic solutions, unlike the case 1T3 — 0, when they only exist as 
a limiting form of the cuspon solutions corresponding to k 0. As far as we know, the model formulated 
in the present work is the first spatially uniform non-degenerate one (i.e., a model with a non- vanishing 
linear part) which yields both cuspons and peakons. 
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1.2 Stability of the solitons and spatiotemporal collapse 

As concerns the dynamical stability of the various solitons in the model (Q) -(^), in this work we limit our- 
selves to direct simulations, as a more rigorous approach, based on numerical analysis of the corresponding 
linear stability-eigenvalue problem [p3| , is technically difficult in the case of cuspons and peakons (results 
of such an analysis, based on the Evans- function technique, will be presented elsewhere). In fact, direct 
simulations of perturbed cuspons and peakons is a hard problem too, but we have concluded that identical 
results concerning the stability are produced (see section 3 below) by high-accuracy finite-difference and 
pseudo-spectral methods (each being implemented in more than one particular form), which lends the 
results credibility. A general conclusion is that the regular solitons are always stable. As for the cuspons 
and peakons, they may be either stable or unstable. 

If the cusp is strong enough, the numerical results presented below demonstrate that the instability 
of a cuspon initiates formation of a genuine singularity, i.e., onset of a spatiotemporal collapse in 
the present one-dimensional model. Before proceeding to the consideration of solitons in the following 
sections, it is relevant to discuss collapse phenomenon in some detail. 

A simple virial-type estimate for the possibility of the collapse can be done, assuming that the field 
focuses itself in a narrow spot with a size (t), amplitude H(i), and a characteristic value K{t) of the field's 
wavenumber |2^. The conservation of the norm (||) imposes a restriction H^L N, i.e., L ~ iV/H^. 
Next, the self-focusing part (||) of the Hamiltonian (||), which drives the collapse, can be estimated as 

i?focus ~ ~ -A^H^. (13) 

On the other hand, the collapse can be checked by the gradient term ( |l0| ) in the full Hamiltonian, that, 
in the same approximation, can be estimated as TJgrad ~ 'i^^KL ^ NK. Further, Eqs. (||) - (^ suggest 
an estimate AT for a characteristic wavenumber of the wave field [the same estimate for K follows 

from an expression ( pT] ) for the exact stationary-soliton solution given below], thus we have -ffgiad ^ A^H^. 
Comparing this with the expression (|l3|) , one concludes that the parts of the Hamiltonian promoting and 
inhibiting the collapse scale the same way as H — + oo (or L ^ 0), hence a weak collapse | |2^ may be 
possible (but does not necessarily take place) in systems of the present type. We stress that, in one- 
dimensional models of GS studied thus far and based on GMTM, collapse has never been reported. The 
real existence of the collapse in the present one-dimensional three-wave GS model, which will be shown 
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in detail below as a result of numerical simulations, is therefore a novel dynamical feature, and it seems 
quite natural that cuspons and peakons, in the case when they are unstable, play the role of catalysts 
stimulating the onset of the collapse. The possibility of a real collapse in a ID system is quite interesting 
by itself, and also because experimental observation of spatio-temporal self-focusing in nonlinear optical 
media is a subject of considerable interest, see, e.g., Ref. [ p5[ . 

2 Analytical solutions 

2.1 The dispersion relation 

The first step in the investigation of the system is to understand its linear spectrum. Substituting 
ui,2,3 ~ exp{ikx — icot) into Eqs. ^ and omitting nonlinear terms, we arrive at a dispersion equation, 

(w^ - fc2 - - ujo) = 2k^{lj ~ 1). (14) 

If K = 0, the third wave decouples, and the coupling between the first two waves produces a commonly 
known gap, so that the solutions to Eq. (|lj) are wi_2 — ±\/l + and = wo- If k 7^ 0, the spectrum 
can be easily understood by treating k as a small parameter. However, the following analysis is valid for 
all values of k in the range < < 1. 

First, consider the situation when fc = 0. Three solutions of Eq. (^4|) arc then 

= \,uj = u± = (0)0 - l)/2 ± ^J{u}Q + 1)2/4 + 2k2. (15) 

It can be easily shown that < minjo^o, —1} < maxjcjo, ~1} < so that one always has cij_ < — 1, 
while LJ+ < 1 if < 1 — k^, and cij+ > 1 if o^g > 1 — ■ Next, it is readily seen that, as k'^ — > cx3, 
either cj^ « fc^, or k, loq. Each branch of the dispersion relation generated by Eq. (^^ is a monotonic 
function of fc^. Generic examples of the spectrum are shown in Fig. 2, where the panels (a) and (b) 
pertain, respectively, to the cases < 1 — with lo^ < 1, and t^o > 1 with > 1. The intermediate 
case, 1 — < < 1, is similar to that shown in panel (a), but with the points lo^ and 1 at fc = 
interchanged. When ujq < 1, the upper gap in the spectrum is min{w+,l} < uj < max{tj+, 1}, while 
the lower gap is lj_ < uj < ujq. When luq > 1, the upper gap is wo < < (^+, and the lower one is 
cj_ < a; < 1. 
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2.2 Gap solitons 

The next step is to search for GS solutions to the full nonlinear system. In this work, we confine ourselves 
to the case of zero- velocity GS, substituting into Eqs. - (||) 

Un{x,t) = Un{x) e'xp{—iu!t) , n — 1,2,3, (16) 

where it is assumed that the soliton's frequency uj belongs to one of the gaps. In fact, even the description 
of zero-velocity solitons is quite complicated. Note, however, that if one sets k = in Eqs. (Q) - (^), 
keeping nonlinear XPM couplings between the first two and third waves, the gap which exists in the two- 
wave GMT model remains unchanged, and the corresponding family of GS solutions does not essentially 
alter, in accord with the principle that nonlinear couplings cannot alter gaps or open a new one if the 
linear coupling is absent [|T^; nevertheless, the situation is essentially different if k is vanishingly small, 
but not exactly equal to zero, see below. 

The substitution of (H) into Eqs. (|) and (|) leads to a system of two ordinary differential equations 
for Ui{x) and U2{x), and an algebraic relation for U^lx), 

iU[ =cjC/i + C/2 + K[/3 + a(aai|t/iP + a|C/2p + |C/3nC/i, (17) 
- iU'2 = UU2 + C/i + Kf/3 + a {aai\U2\^ + a\Ui\^ + IC/3P) U2, (18) 
(c^o - ^ + fT3|C/3p + a|C/i|2 + a\U2\'')U^ = k(C/i + U2), (19) 

where the prime stands for d/dx. To solve these equations, we substitute Ui^2 = ^1,2(2:) exp (z(/)i_2(2:)) 
with real An and After substituting the expression ( p^ ) into equations ( p^ , p^ , and some simple 
manipulations, it can be found that {^A\ — A2)' = and {(pi + (j)2)' = 0. Using the condition that the 
soliton fields vanish at infinity, we immediately conclude that 

Al{x) ^ Alix) = S{x); (20) 

as for the constant value of 0i +02, it may be set equal to zero without loss of generality, so that 
(l)i{x) = —(j)2{x) = 0(x)/2, where (j){x) is the relative phase of the two fields. After this, we obtain two 
equations for S{x) and (/)(x) from Eqs. ([l7| ) and (^, 

(j)' = -2LU-2cos(l3~2a^{l + ai)S-S-^Ui{ujQ-uj-a3Ui) , (21) 
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S" = -25sin0-2KV5t/3sin(0/2), (22) 
and Eq. (|9|) for the third wave U3 takes the form of a cubic algebraic equation 

Us {ujQ~uj~ 2aS ~ CTslC^sH = 2kv^cos (0/2) , (23) 

from which it foUows that t/3 is a real-valued function. 

This analytical consideration can be readily extended for more general equations (^) and (||) that do 
not assume the symmetry between the waves ui and U2, i.e., with the group- velocity terms in the equations 
altered as in Eq. (^. In particular, the relation ( po|) is then replaced by ciA^(x) — £2^2(2;) = S{x). 
The subsequent analysis is similar to that above, and leads to results for the asymmetric model that are 
qualitatively similar to those presented below for the symmetric case. 

Equations (|2l| ) and (|2^) have a Hamiltonian structure, as they can be represented in the form 

dS _ dH d(j) _ dH ^2^^ 
dx dcj) ' dx dS ' 

with the Hamiltonian 

H = 2Scos(l) + c? (1 + CTi) + 2ujS + C/| [ujq-^- 2aS) - ^^3^/3 , (25) 

which is precisely a reduction of the Hamiltonian (||) of the original system (|^) - (|^) for the solutions of 
the present type. Note that H is here regarded as a function of S and (j), and the relation ( p3|) is regarded 
as determining C/3 in terms of S and </>. We stress that the dependence U^^SjCf)) was taken into account 
when deriving the Hamiltonian representation (|2^). 

For soliton solutions, the boundary conditions at a; = ±00 yield H ~ so that the solutions can be 
obtained in an implicit form, 

2Scos(j) + a'^{l + ai)S'^ + 2LoS+Ui{LOa-uj-2aS)-{3/2)a3Ul = 0. (26) 



In principle, one can use the relations (|2^) and (g6|) to eliminate U3 and (/> and so obtain a single equation 
for S. However, this is not easily done unless as — Q (no SPM term in Eq. (ref3)), and so we proceed to 
examine this special, but important, case first. Note that the no-SPM case also plays an important role 
for GMTM, which is exactly integrable by means of the inverse scattering transform just in this case Q]. 
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2.3 Cuspons, the case = 

Setting (73 = makes it possible to solve Eq. (|2^) for 3 explicitly in terms of S and (f>, 



2kVScos{cI)/2) 
LjQ — uj — Zab 

For simplicity, we also set cti = in Eqs. (|^) and (||) and subsequent equations, although the latter 
assumption is not crucially important for the analysis developed below. Indeed, the analysis is based on 
the fact that the field U3 can be explicitly eliminated by means of Eq. (|2^), which is not affected by cti. 
If (Ti is kept in the system, it merely renormalizes some coefficients in the formulas derived below. 

At the next step, one can also eliminate (p, using Eqs. ) and (^), to derive a single equation for 

S, 

[dS/dxf = AS'^F{S), (28) 

1 



F{S)^{l-u;-^a'S) 



2(1 + ^— -) - {l-^-^\a^S) 



(29) 



The function F{S) has cither one or three real zeros S^. One is 

Soi = 2 (l-c^)/a^ (30) 
and the remaining two, if they exist, are real roots of the quadratic equation, 

(2 + 2w + a^So){uJa ~uj- 2aSo) + Ak^ = 0. (31) 



Only the smallest positive real root of Eq. (|31|), to be denoted S02 (if such exists), will be relevant below. 
Note, incidentally, that F{S) cannot have double roots. It is easy to see that a consequence of this fact is 
that Eq. ( p^ ) cannot generate kink solutions, which have different limits as a; — > ±00, for both of which 
the right-hand-side of (|8|) must have a double zero. 

For a bright-soliton solution of (p8|), we need first that F(Q) > (in this work, we do not consider 
dark solitons, nor "anti-dark" solitons, i.e., solitons on top of a finite-amplitude flat background, a reason 
being that there is little chance that the flat background would be modulationally stable). Comparing 
the condition F{0) > with the expressions given in Section 2.1 for the gaps in the linear spectrum, 
it is readily shown that this condition is exactly equivalent to requiring that lo belongs to either the 
upper or the lower gap of the linear spectrum. We note that the coupling to the third wave gives rise to 



nonlinearity of the rational type in the expression (29), despite the fact that the underlying system 
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contains only cubic polynomial nonlinear terms. Even if the coupling constant k is small, it is clear 
that the rational nonlinearity may produce a strong effect in a vicinity of a critical value of the squared 
amplitude at which the denominator in the expression ( p9| ) vanishes, 

5„ = (c^o-c^)/2a, (32) 

where one must have a{ujQ — lj) > (otherwise, this critical value is not relevant). 

If S'cr > 0, the structure of the soliton crucially depends on whether, with an increase of 5, the 
function F{S) defined by Eq. (|2^ ) first reaches zero at 5 = 5*0 > (i.e. either S — Sm or S = Sq2, 
whichever is the smaller positive value), or, instead, it first reaches the singularity at 5 = Scr, i.e., whether 
< So < Sci, or < S'cr < Sq. In the former case, the existence of Sa plays no role, and the soliton 
is a regular one, having the amplitude ^/Sq. This soliton may be regarded as obtained by a smooth 
deformation from the usual GS known in GMTM at k = 0. 

In the case < S'cr < So, as the soliton cannot have an amplitude larger than \/S^, the amplitude 
takes this critical value. The soliton is singular in this case, being a cuspon (see details below), but, 
nevertheless, it is an absolutely relevant solution. The remaining possibilities are that either Scr < 
and So > 0, or vice versa; then the soliton may only be, respectively, regular or singular. Of course no 
soliton exists if both So and Scr are negative. Further, using the symmetries of the equations, it is readily 
shown that for all these soliton solutions, S{x) is symmetric about its center, which may be set at x — 0, 
that is, S{x) is an even function of x. For the cuspon solutions, and for those regular solutions whose 
squared amplitude is Soi, it can also be shown that the phase variable ip{x) = 4>{x) — n and Uslx) are 
odd functions of x, while for those regular solutions whose squared amplitude is S02 the phase variable 
(f>{x) and U'i{x) are, respectively, odd and even functions of x. 

It is now necessary to determine which parameter combinations in the set (w, wo, a) permit the options 
described above. The most interesting case occurs when loq > uj (so that lj belongs to the lower gap, 
see Fig. 2) and a > (the latter condition always holds in the applications to nonlinear optics). In this 
case, it can be shown that the root S02 of Eq. ( |3l| ) is not relevant, and the options are determined by 
the competition between S'oi and Scr- The soliton is a cuspon (0 < Scr < Soi) if 

a(wo-w) < 4(1 -w). (33) 

In effect, the condition (|3^ ) sets an upper bound on a for given loq and lo. In particular, this condition 
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is always satisfied if < a < 4. 

If, on the other hand, the condition ( ^3| ) does not hold (i.e., < Sqi < S'cr), we obtain a regular 
soliton. In a less physically relevant case, when again loq > lu but a < 0, cuspons do not occur (as this 
time 5cr < 0, see Eq. (^, and only regular solitons may exist. 

Next we proceed to the case ujq < oj, so that uj is located in the upper gap of the linear spectrum. For 
a > 0, we have S'cr < 0, hence only regular solitons may occur, and indeed in this case there is always 
at least one positive root 5*0, so a regular soliton does exist. If a < 0, then we have S'cr > 0, but if 
llIq < 1 — k'^ (when also cu < 1), there is at least one positive root So < Scr; thus, only a regular soliton 
can exist in this case too. On the other hand, if a < and uq > 1 — (and then a; > 1), there are no 
positive roots Sq, and so only cuspons occur. 

Let us now turn to a detailed description of the cuspon's local structure near its center, when S is 
close to Scr. From the above analysis, one sees that cuspons occur whenever u; lies in the lower gap, 
with ujQ > uj and a > 0, so that the criterion ( ^3| ) is satisfied, or when lu lies in the upper gap with 
1 — < ujQ < Lo and a < 0. To analyze the structure of the cuspon, we first note that, as it follows from 
Eq. ([2^), one has cos0 = —1 (i.e., = tt) when S — Scr, which suggest to set 

S„-S = 5 - k'^R, 1 + cos(I) = S ■ p, (34) 

where ^ is a small positive parameter, and the stretched variables R and p are positive. At the leading 
order in 6, it then follows from Eq. ( p6|) that p = poR, where 

po = a^(Soi - S'cr). (35) 

As it follows from the above analysis, po is always positive for a cuspon. We also stretch the spatial 
coordinate, defining x = S'^^^K.^y, the soliton center being at x = 0. Since S{x) is an even function of x, 
it is sufficient to set a; > in this analysis. Then, on substituting the first relation from Eq. ( |3^ ) into 
Eq. ([2^), we get, to the leading order in S, an equation 

R{dR/dyf = PoSl/a^ = K^, (36) 

so that 

R = {?,Ky/2f'^ . (37) 
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Note that in the original unstretched variables, the relation (|3^ ) shows that, near the cusp, 

5cr - S{x) w {3KKx/2f/^, (38) 

dS/dx « (2/3)1/''' (if/^)'/^ • (39) 
and it follows from Eq. ( p7| ) that C/3 is unbounded near the cusp, 

U3 « iScr/a){2apoKy3,ixy/\ (40) 



In particular, Eq. (39) implies that, as Kk decreases, the cusp gets localized in a narrow region where 
^K^K^ (outside this region, \dS/dx\ is bounded and shows no cusp). Note that this limit can be 
obtained either as 0, or as po ^ (recall that po is defined in Eq. (|3^)). 

An example of the cuspon is shown in Fig. 3. Although the first derivative in the cuspon is singular 
at its center, as it follows from Eq. ( ^9[ ) [see also Fig. 3(a)], and its U3 component diverges at x 
as per Eq. (^), it is easily verified that the value of the Hamiltonian (||) [and, obviously, the value of 
the norm (H) too] is finite for the cuspon solution. These solitons are similar to cuspons found as exact 
solutions to the Camassa-Holm (CH) equation [jl^, which have a singularity of the type \x\^^'^ or 
|a;|^/^ as |a;| —* 0, cf. Eqs. (|3|) and (|39|). The CH equation is integrable, and it is degenerate in the 
sense that it has no linear terms except for du/dt (which makes the existence of the solution with a cusp 
singularity possible). Our three- wave system (|^) - (H) is not degenerate in that sense; nevertheless, the 
cuspon solitons are possible in it because of the model's multicomponent structure: the elimination of 
the third component generates the non-polynomial nonlincarity in Eqs. (|l7|), ([l8|), and, finally, in Eqs. 
(|2^ ) and (|2^), which gives rise to the cusp. It is noteworthy that, as well as the CH model, ours gives 
rise to two different coexisting families of solitons, viz., regular ones and cuspons. It will be shown below 
that the solitons of both types may be stable. 

Of course, the presence of the singularities in U3{x) and dS/dx at x ^ suggests that higher-order 
terms, such as the higher-order dispersion, should be taken into regard in this case. The fact that the 
cuspon's Hamiltonian converges despite these singularities, as well as a direct analysis, suggest that such 
higher-order terms will smooth the shape of the cuspon in a very narrow layer for small x, allowing for 
large but not diverging values of the fields. However, the small higher-order terms will not essentially 
alter the global shape of the cuspons. In the next section we will show that, in fact, the genuine generic 
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singular solitons are (in the presence of the SPM terms) peakons, for which the singularities are much 
weaker, hence the latter issue is still less significant. Besides that, it appears to be an issue of principal 
interest to understand what types of solitons the system may generate without intrinsic dispersion (cf. 
the situation for the traditional GMTM, in which the spectrum of soliton solutions is completely altered 
by the addition of intrinsic dispersion p^). 

In the special case k ^ 1, when the third component is weakly coupled to the first two ones in 
the linear approximation (in terms of the optical model represented by Fig. 1, it is the case when the 
subgratings shown by the dashed lines are very weak), straightforward inspection of the above results 
shows that the cuspons look like peakons; that is, except for the above-mentioned narrow region of the 
width \x\ ^ K^, where the cusp is located, they have the shape of a soliton with a discontinuity in the first 
derivative of S{x) and a jump in the phase 4>{x), which are the defining features of peakons (||l^, |20|). 
A principal difference of true peakons from cuspons is that the first derivative does not diverge inside a 
peakon, but is of course, discontinuous. 

An important result of our analysis is that the family of solitons obtained in the limit k — s- is 
drastically different from that in the model where one sets k = from the very beginning. In particular, 
in the most relevant case, with loq > lu and a > 0, the family corresponding to k — > contains regular 
solitons whose amplitude is smaller than \/Sci-; however, the solitons whose amplitude at k = is larger 
than-v/S'cr, i.e., the ones whose frequencies belong to the range (|3^ ) [note that the definition of Sa does 
not depend on k at all, see Eq. (|3^)], are replaced by the peakons which are constructed in a very simple 
way: drop the part of the usual soliton above the critical level S — Scr, and bring together the two 
symmetric parts which remain below the critical level, see Fig. 3(b). 

It is interesting that peakons are known as exact solutions to a version of the integrable CH equation 
slightly different from that which gives rise to the cuspons. As well as in the present system, in that 
equation the peakons coexist with regular solitons [ pO[ . In the next subsection, we demonstrate that 
the peakons, which are found only as limit-form solutions in the no-SPM case CT3 = 0, become generic 
solutions in the case 0-3 ^ 0. 
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2.4 Peakons, the case as 7^ 

A natural question is whether the cuspon solutions are structurally stable, i.e., if they will persist on 
inclusion of terms that were absent in the analysis presented above (the other type of the stability, viz., 
dynamical stability against small initial perturbations, will be considered in the next section). Here, we 
address this issue by restoring the SPM term in Eq. (|^), that is, we now set 0-3 7^ 0, but assume that it 
is a small parameter. Note that, in the application to nonlinear optics, one should expect that 173 > 0, 
but there is no such a restriction on the sign of (T3 in the application to the flow of a density-stratified 
fluid. We still keep ai = 0, as the inclusion of the corresponding SPM terms in Eqs. (|^) and (^) amounts 
to trivial changes both in the above analysis, and in that presented below. On the other hand, we show 
below that the inclusion of the SPM term in Eq. (||) is a structural perturbation which drastically changes 
the character of the soliton solutions. 

In view of the above results concerning the cuspons, we restrict our discussion here to the most 
interesting case when S{x) is an even function of x, while tpi^) — (f'i^) ~ ^ ^-nd U3{x) are odd functions. 
In principle, one can use the relations (|2^) and ( p6|) to eliminate (j) and U3 and so obtain a single equation 
for S (a counterpart to Eq. (p8|)), as it was done above when = 0. However, when (73 7^ 0, it is not 
possible to do this explicitly. Instead, we shall develop an asymptotic analysis valid for a; — > 0, which 



will be combined with results obtained by direct numerical integration of Eqs. (21) and (|22|), subject of 
course to the constraints ( |2^ ) and (p6|). Since singularities only arise at the center of the soliton (i.e., 
at a; = 0) when (T3 = 0, it is clear that the introduction of a small (T3 7^ will produce only a small 
deformation of the soliton solution in the region where x is bounded away from zero. 

First, we consider regular solitons. Because the left-hand side of Eq. ( p3| ) is not singular at any x, 
including the point a; = 0, when (T3 — 0, we expect that regular solitons survive a perturbation induced 
by fs 7^ 0. Indeed, if there exists a regular soliton, with Sq = S{x = 0), and 4>{x = 0) = vr and 
U^ix = 0) = 0, it follows from Eq. ( ^6|) that the soliton's amplitude remains exactly the same as it was 
for (T3 = 0, due to the fact that the regular soliton has {/^{x = 0) = 0. 

Next, we turn to the possibility of singular solutions, that is, cuspons or peakons. Since we are 
assuming that Sq — S{x — 0) is finite, and that 4>{x ~ Q) ~ n, it immediately follows from Eq. ([2^) that 



when CT3 7^ 0, U3 must remain finite for all x, taking some value Uq =/= 0, say, as x — > +0. Since U3 is an 
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odd function of x, and Uq ^ 0, there must be a discontinuity in U3 at x — 0, i.e., a jump from Uq to —Uq. 
This feature is in marked contrast to the cuspons for which C/3 is infinite at the center, see Eq. (|40|). 
Further, it then follows from Eq. ( p2| ) that, as a; — > 0, there is also a discontinuity in dS/dx, with a jump 
from 2k[/o-\/^ to —2kUo^/So- Hence, if we can find soliton solutions of this type, with Uq ^ 0, they 
are necessarily peakons, and we infer that cuspons do not survive the structural perturbation induced by 

(73^0. 

Further, if we assume that Uq 7^ 0, then Eq. (|2^), taken in the limit x ^ 0, immediately shows that 

2a(5'cr - 5*0) = asU^ (41) 



(recall that Sct is defined by Eq. (p2|)). Next, the Hamiltonian relation (26), also taken in the limit 
X ^ 0, shows that 

- ^So - a^SoiSc, - So) = lasU^, (42) 



where we have used Eq. (|4l|) [recall that po is defined by Eq. (|35|)]. Elimination of from Eqs. ( |4l| ) and 
( p^ ) yields a quadratic equation for Sq , whose positive roots represent the possible values of the peakon's 
amplitude. 

We recall that for a cuspon which exists at (T3 = one has po > 0, i.e., the amplitude of the 
corresponding formal regular soliton exceeds the critical value of the amplitude, see Eq. (p5|). Then, if 
we retain the condition po > 0, it immediately follows from Eqs. (^) and (^2|) that no peakons may exist 
if the SPM coefficient in Eq. (||) is positive, 0-3 > 0. Indeed, Eq. (|4|) shows that 5cr — S'o > if (T3 > 0, 
which, along with po > 0, leads to a contradiction in the relation (^). 

Further, it is easy to see that a general condition for the existence of peakons following from Eqs. 
(|l|) and dH) is 

<y?,PQ < , (43) 

hence peakons are possible if 173 < 0, or if we keep 173 > but allow po < 0. In the remainder of this 
section, we will show that peakons may exist only if po > 0. Hence, it follows from the necessary condition 
( ^ ) that peakons may indeed be possible solely in the case (T3 < 0. On the other hand, regular solitons 
do exist in the case 173 > (i.e., in particular, in nonlinear-optics models), as they have Uq — 0, hence 
neither Eq. ([4l| ) nor its consequence in the form of the inequality (Bs) apply to regular solitons. The 
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existence of (stable) peakons for 0-3 < 0, and of (also stable) regular solitons for > will be confirmed 
by direct numerical results presented in the next section. 

To obtain a necessary condition (which will take the form of po > 0) for the existence of the peakons, 
we notice that the existence of any solitary wave implies the presence of closed dynamical trajectories 
in the phase plane of the corresponding dynamical system, which is based on the ordinary differential 
equations (|2^) and (p^), supplemented by the constraint (p3|). Further, at least one stable fixed point 
(FP) must exist inside such closed trajectories, therefore the existence of such a stable FP is a necessary 
condition for the existence of any solitary wave. 

The FPs are found by equating to zero the right-hand sides of Eq. (|2^ ) and (^), which together 



with Eq. (23) give three equations for the three coordinates ^, S and C/3 of the FP. First of all, one can 
find a trivial unstable FP of the dynamical system, 

UJ + / {UJO - Uj) 

cos0 = --— — Y77 7' ^^0, 

which does not depend on (T3. Then, three nontrivial FPs can be found, with their coordinates i^*, 5* 
and C/3* given by the following expressions: 

4^^=-, Si'^ ='-f^lsou U^'J=0, (44) 
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TT, (2-^3)5^ =2S'„.- y^oi, (2-(T3) aC/if =po~a''S,,, (45) 
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(2 - a3)5f ^ = 2S„ - ^a;Soi +'—, (2 - a^) [aU^l' \ = po - a'S„ - a^K' , 

cos (4'V2) = /V^ ' (46) 

where the superscript is a number label for the FP. To be specific, we now consider the case of most 
interest, when both 5*01 > and Sci- > . In this case, the FP given by Eqs. (Q ) exists for all CT3 and 
all pq. However, for small 0-3 (in fact 0-3 < 2 is enough) and small k, the FPs given by Eqs. (||) and (||) 
exist only when po > 0. Indeed, they exist only for po > a^Soi and po > a'^S'oi + k^, respectively, or, on 
using the definition (|5| ) of po, when 6*01 > 2S'cr and 5oi > 25cr + K^/a, respectively. 

Let us first suppose that po < 0. Then there is only the single non-trivial FP, namely the one given 
by Eqs. (p4[). This FP is clearly associated with the regular solitons, whose amplitude at the crest is 5oi. 
Hence, we infer that for pQ < there are no other solitary-wave solutions, and in particular, no peakons 
(and no cuspons when = either, in accordance with what we have already found in subsection 2.3 
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above). When combined with the necessary condition (^) for the existence of peakons, we infer that 
there are no peakons when (T3 > 0, thus excluding peakons from apphcations to the nonhnear-optics 
models, where this SPM coefhcient is positive. However, peakons may occur in density-stratified fluid 
flows, where there is no inherent restriction on the sign of (73 . This case is considered below, but first we 
note that in the case po < and CT3 > (which includes the applications to nonlinear optics), the same 
arguments suggest that there may be periodic solutions with a peakon-type discontinuity at the crests; 
indeed, our numerical solutions of the system (^ , ^ ) show that this is the case. 

Next, we suppose that po > 0. First, if 5*01 < 2Sci then there is again the single non-trivial FP given 
by (^4|). But now, by analogy with the existence of cuspons when po > and — , , we infer that the 
solitary wave solution which is associated with this fixed point is a peakon, whose squared amplitude Sq 
for small (73 is close to Scr, while the FP has si^^ — S'oi/2 < S'cr- 



If, on the other hand, 5oi > 2S'cr, the FPs given by Eqs.([45D and (46) become available as well. We 
now infer that the peakon solitary-wave solution continues to exist, and for sufhciently small (T3 and k 
it is associated with the FP given by Eq. (45). Although Eq. (Esl) implies that S'i^'' w Sa, and the 



peakon's squared amplitude Sq, determined by Eqs. ( p| ) and (^2|), is also approximately equal to Scr, 
we nevertheless have 6*0 > S'*^'' as required. Note that, in the present case, the FPs given by Eqs. (^) 
and ( p6| ) lie outside the peakon's homoclinic orbit. In Fig. 4, we show a plot of a typical peakon obtained, 
in this case, by numerical solution Eqs. (|2l|) and (p2[). 



3 Numerical results 
3.1 Simulation techniques 

The objectives of direct numerical simulations of the underlying equations (^ - (^ were to check the 
dynamical stability of regular solitons, cuspons, and peakons in the case = 0, and the existence and 
stability of peakons in the more general case, 173 7^ 0. Both finite-difference and pseudo-spectral numerical 
methods have been used, in order to check that the same results are obtained by methods of both types. 
We used semi-implicit Crank-Nicholson schemes, in which the nonlinear terms were treated by means of 
the Adams-Bashforth method. 

The presence of singularities required a careful treatment of cuspon and peakon solutions. To avoid 
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numerical instabilities due to discontinuities, we found it, sometimes, beneficial to add an artificial weak 
high-wavenumber viscosity to the pseudospectral code. This was done by adding linear damping terms to 
the right-hand side of Eqs. (^, (^) and (|6|), which have the form — ii/(fc)fc^u„ in the Fourier representation, 
where u„ is the Fourier transform of w„ {n ~ 1, 2, 3). The high-pass filter viscosity i^(fc) suppresses only 
high wavenumbers and does not act on others. In particular, we chose 

0, if|fc|<^if 
v(k)^\ ^(1M_5), \i^^K<\k\<lK , 
7?, if > \K 

where K is the largest wavenumber in the actual numerical scheme, and 77 is a small viscosity coefficient. 
We have found that 77 ~ 10^^ was sufficient to avoid Gibbs' phenomenon in long-time simulations. 

When instabilities occur at a singular point (cusp or peak), it is hard to determine whether the 
instability is a real one, or a numerical artifact. Therefore, we checked the results by means of a finite- 
difference code which used an adaptive staggered grid; motivated by the analysis of the vicinity of the 
point a: = reported above, we introduced the variable ^ = x^l'^ to define an adaptive grid, and also 
redefined C/3 = \/\lJ-i. In these variables, the cusp becomes a regular point. This approach was solely 
used to check the possible occurrence of numerical instabilities. 

In the following subsections we present typical examples of the numerical results for both cases con- 
sidered above, viz., (73 = and (T3 < 0, when, respectively, cuspons and peakons are expected. 

3.2 The case (T3 = 

First, we report results obtained for the stability of regular solitary waves in the case (T3 = 0. As initial 
configurations, we used the corresponding stationary solutions to Eqs. and (|2^). To test the stability 
of the regular solitary waves, we added small perturbations to them. As could be anticipated, the regular 
solitary wave sheds off a small dispersive wave train and relaxes to a stationary soliton, see Fig. 5 (for a 
more detailed illustration of the generation of small radiated waves by a soliton, see also Fig. 8 below). 
If, however, a regular soliton is taken as an initial condition for parameter values inside, but close to the 
border of the cuspon region, it does not become unstable in this slightly modified section of parameter 
space (which only supports cuspons), but instead this soliton exhibits persistent internal vibrations, see 
an example in Fig. 6. These and many other simulations clearly show that the regular soliton is always 
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stable, and, close to the parameter border with cuspons, it has a persistent internal mode. 

It was shown analytically above that Eqs. (|l^) and ( [l8| ) ( with 0-3 = 0) support peakons when po > 
and poK^ is very small. Direct simulations show that peakons do exist in this case, and they may be 
either unstable or stable. In the case when they are unstable, a high-wavenumber instability develops 
around the central peak. In Fig. 7, we display the time evolution of a typical stable peakon. 

Next, we look at what happens if we take a regular soliton as an initial condition in a section of the 
parameter space which supports only stable peakons. This enables us to study the competition of the 
structural stability of regular solitons (as confirmed in Fig. 6) and the stability of peakons. The initial 
condition is taken as a stationary regular soliton in the parameter region (close to the boundary of the 
peakon region) with po < 0, whereas the simulations are run for values of the parameters corresponding 
to Po > 0, which only admits peakons and excludes regular solitons. Unlike the case shown in Fig. 6, the 
time evolution now does not exhibit internal vibrations. Instead, the pulse slowly decays into radiation. 
This outcome can be explained by the fact that the peakon's norm (see Eq.(||)) turns out to be larger 
than that of the initial pulse in this case, hence its rearrangement into a stable peakon is not possible. 
An essential result revealed by the simulations is that cuspons may also be stable, a typical example 
being displayed in Fig. 8. In this figure, one can see a small shock wave which is initially generated 
at the cuspon's crest. It seems plausible that this shock wave is generated by some initial perturbation 
which could be a result of the finite mesh size in the finite-difference numerical scheme employed for the 
simulations. In fact, the emission of a small-amplitude shock wave is quite a typical way of the relaxation 
of both cuspons and peakons to a final stable state. To make it sure that the shock wave is not an artifact 
generated by the numerical scheme, we have checked that its shape does not change with the increase of 
the numerical accuracy. 

To further test the stability of the cuspons and peakons, in many cases we allowed the initially 
generated shock wave to re-enter the integration domain (due to periodic boundary conditions used in 
the simulations) and interact again with the cuspon or peakon. As a result, the stability of the solitons of 
these types has been additionally confirmed. An example of the spatial profile of the cuspon established 
after a long evolution is shown in Fig. 9. Both the stability of the cuspon, and the presence of a tiny 
shock wave are evident in the figure. 

However, unlike the regular solitons, which were found to be always stable, the cuspons are sometimes 
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unstable. Typically, their instability triggers the onset of spatiotemporal collapse, i.e., formation of a 
singularity in a finite time (see a discussion of the feasible collapse in systems of the present type, given 
in the Introduction). Simulations of the collapse were possible with the use of an adaptive grid. A typical 
example of the collapse is shown in Fig. 10, where the inset shows that (within the numerical accuracy 
available) the amplitude of the collapsing pulse indeed diverges in a finite time. 

However, the collapse is not the only possible outcome of the instability. In some other cases, which 
are not displayed here, the instability of peakons could be quite weak, giving rise to their rearrangement 
into regular solitons by the shedding of a small amount of radiation. 

3.3 The case (73 7^ 

The predictions of the analysis developed above for the most general case, when the SPM terms are 
present in the model (0-3 ^ 0), were also checked against direct simulations. As a result, we have found, 
in accord with the predictions, that only regular solitons exist in the case 0-3 > 0, while in the case as < 
both regular solitons and peakons have been found as generic solutions. Further simulations, details of 
which are not shown here, demonstrate that both regular solitons and peakons are stable in this case. 

4 Conclusion 

In this work, we have introduced a generic model of three waves coupled by linear and nonlinear terms, 
which describes a situation when three dispersion curves are close to an intersection at one point. The 
model was cast into the form of a system of two waves with opposite group velocities that, by itself, gives 
rise to the usual gap solitons, which is further coupled to a third wave with zero group velocity (in the 
laboratory reference frame). Situations of this type are quite generic, being realizable in various models of 
nonlinear optics, density-stratified fluid flows, and in other physical contexts. In particular, two versions 
(temporal and spatial) of a nonlinear-optical model, which is based on a waveguide carrying the triple 
spatial Bragg grating, have been elaborated in the Introduction. Our consideration was focussed on zero- 
velocity solitons. In a special case when the self-phase modulation (SPM) is absent in the equation for 
the third wave, soliton solutions were found in an exact form. It was shown that there are two coexisting 
generic families of solitons: regular solitons and cuspons. In the special case when the coefficient of the 
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linear coupling between the first two waves and the third one vanishes, cuspons are replaced by peakons. 
Direct simulations have demonstrated that the regular solitons are stable (in the case when the regular 
soliton is close to the border of the cuspon region, it has a persistent internal mode). The cuspons and 
peakons may be both stable and unstable. If they are unstable, they either shed off some radiation and 
rearrange themselves into regular solitons, or, in most typical cases, the development of the cuspon's 
instability initiates onset of spatiotemporal collapse. Actually, the present system gives the first explicit 
example of collapse in one-dimensional gap-soliton models. 

The most general version of the model, which includes the SPM term in the equation for the third 
wave, has also been considered. Analysis shows that cuspons cannot exist in this case, i.e., cuspons, 
although being possibly dynamically stable, are structurally unstable. However, depending on the signs 
of the SPM coefficient, and some combination of the system's parameters, it was shown that a generic 
family of peakon solutions may exist instead. In accord with this prediction, the peakons have been found 
in direct simulations. The peakons, as well as the regular solitons, are stable in the system including the 
SPM term. We stress that peakons are physical solutions, as they have all their field components and 
their first derivatives finite. 

The next step in the study of this system should be consideration of moving solitons, which is suggested 
by the well-known fact that the usual two- wave model gives rise to moving gap solitons too j^. However, 
in contrast to the two-wave system, one may expect a drastic difference between the zero-velocity and 
moving solitons in the present three-wave model. This is due to the reappearance of a derivative term 
in Eq. (^), when it is written for a moving soliton, hence solitons which assume a singularity or jump 
in the U3 component, i.e., both cuspons and peakons, cannot then exist. Nevertheless, one may expect 
that slowly moving solitons will have approximately the same form as the cuspons and peakons, with the 
singularity at the central point replaced by a narrow transient layer with a large gradient of the U3 field. 
Detailed analysis of the moving solitons is, however, beyond the scope of this work. 

Acknowledgements 

We would like to thank Tom Bridges, Gianne Derks and Sebastian Reich for valuable discussions. B.A.M. 
and G.A.G. appreciate hospitality of the University of Loughborough (UK). The work of G.A.G. is 



25 



supported by a European Commission Grant, contract number HPRN-CT-2000-00113, for the Research 
Training Network Mechanics and Symmetry in Europe (MASIE) 



26 



References 

[1] CM. de Sterke and J.E. Sipe, in Progress in Optics, edited by E. Wolf (Elsevier, North-Holland, 
1994), Vol. XXXIII, Chap. 3, p. 205; A.B. Aceves, Chaos 10, 584 (2000). 

[2] G. Kurizki, A.E. Kozhekin, T. Opatrny, and B.A. Malomed, in Progress in Optics, edited by E. Wolf 
(Elsevier, North-Holland, 2001), Vol. 42, Chap. 2, p. 93. 

[3] R. Grimshaw, J. He, and B.A. Malomed, Physica D 113, 26 (1998). 

[4] D.N. Ghristodoulides and R.I. Joseph, Phys. Rev. Lett. 62, 1746 (1989); A. Aceves and S. Wabnitz, 
Phys. Lett. A 141, 37 (1989). 

[5] N.G.R. Broderick, D.J. Richardson, R.I. Laming, and M. Ibsen, Opt. Lett. 23, 328 (1998). 

[6] B.J. Eggleton, R.E. Slusher, CM. de Sterke, P.A. Krug, and J.E. Sipe, Phys. Rev. Lett. 76, 1627 
(1996); CM. de Sterke, B.J. Eggleton, and P.A. Krug, J. Lightwave Technol. 15, 1494 (1997). 

[7] R. Grimshaw and B.A. Malomed, Phys. Rev. Lett. 72, 949 (1994). 

[8] J. Gear and R. Grimshaw, R. Stud. Appl. Math. 70, 235 (1984). 

[9] B.J. Eggleton BJ, A.K. Ahuja, K.S. Feder, C. Headley, C. Kerbage, M.D. Mermelstein, J.A. Rogers, 
P. Steinvurzel, P.S. Westbrook, and R.S. Windeler, IEEE J. Selected Topics Quant. Electron. 7, 409 
(2001); see also a focus issue: Photonics Crystal Fiber, in Opt. Express 9, 675 (2001). 

[10] D.G. Salinas, CM. de Sterke, and J.E. Sipe, Opt. Commun. Ill, 105 (1994); CM. de Sterke, D.G. 
Sahnas, and J.E. Sipe, Phys. Rev. E 54, 1969 (1996). 

[11] 7 (3): 409-424 MAY-JUN 2001 

[12] W.C.K. Mak, B.A. Malomed, and P.L. Chu, Phys. Rev. E 58 6708 (1998). 
[13] G.P. Agrawal, Nonlinear Fiber Optics (Academic Press: Boston, 1995). 
[14] R. Grimshaw and B.A. Malomed, Phys. Lett. A 198, 205 (1995). 
[15] H. Aratyn, J. Phys. A 14, 1313 (1981). 



27 



[16] T. Peschel, U. Peschel, F. Lederer, and B.A. Malomed, Phys. Rev. E 55, 4730 (1997); J. Yang, B.A. 
Malomed, and D.J. Kaup, Phys. Rev. Lett. 83, 1958 (1999). 

[17] R. Camassa and D.D. Holm, Phys. Rev. Lett. 71, 1661 (1993). 

[18] R.A. Kraenkel and A. Zenchuk, J. Phys. A 32, 4733 (1999); M.C. Ferreira, R.A. Kraenkel, and A.I. 
Zenchuk, J. Phys. A 32, 8665 (1999). 

[19] A.S. Fokas and B. Fuchssteiner, Physica D 4, 47 (1981). 

[20] R. Beals and D.H. Sattinger, Inverse Problems 15, LI (1999). 

[21] T.Qian and M.Tang, Chaos, Solitons and Fractals 12, 1347 (2001). 

[22] N.G.R. Broderick and CM. de Sterke, Phys. Rev. E 58, 7941 (1998). 

[23] I.V. Barashenkov, D.E. PeHnovsky, and E.V. Zemlyanaya, Phys. Rev. Lett. 80, 5117 (1998); A. De 
Rossi, C. Conti, and S. Trillo, ibid. 81, 85 (1998). 

[24] L. Berge, Phys. Rep. 303, 259 (1998). 

[25] H.S. Eisenberg, R. Morandotti, Y. Silberberg, S. Bar-Ad, D. Ross, and J.S. Aitchison, Phys. Rev. 
Lett. 87, 3902 (2001). 

[26] A.R. Champneys, B. A. Malomed, and M.J. Friedman, Phys. Rev. Lett. 80, 4169 (1998). 



28 



Fig. 1. A schematic representation of the optical model that gives rise to the three linearly coupled 
waves in a x^^^ waveguide. The figure shows either the planar waveguide with the triple Bragg grating, 
in the case of the temporal evolution of the fields, or the transverse cross section of the bulk medium 
of the layered (photonic-crystal-fiber) type, in the case of the spatial evolution along the coordinate z 
perpendicular to the plane of the figure. The triangles formed by bold arrows illustrate how the linear 
couplings between the three waves, whose Poynting vectors are represented by the arrows, are induced 
by the Bragg reflections on the three gratings. The difference between the gratings represented by the 
continuous and dashed lines is in their strength (refractive-index contrast). 

Fig. 2. Dispersion curves produced by Eq. (|l4| ) in the case k = 0.5: (a) t^^o < 1 ^ (b) ojq > 1- The 
dashed line in each panel is lu ^ ojq. The case with 1 — < wq < 1 is similar to the case (a) but with 
the points uJ-^- and 1 at fc = interchanged. 

Fig. 3. The shape of the cuspon for a = 2.0, ujq = 0.1, co = —0.5, and (a) k = 0.5, i.e., in the general 
case, and (b) k = 0.1, i.e., for small k. In the case (b) we also show the usual gap soliton (by the dashed 
line), the part of which above the critical value S — Scr (shown by the dotted line) should be removed 
and the remaining parts brought together to form the peakon corresponding to pqk^ ~* 0. 

Fig. 4. The shape of the peakon for the case when 0-3 < 0, where we plot |J7ip. The parameters are 
173 = —0.01, K = 0.1, a = 2.0, Wo = 0.1, and co = —0.5. In this case, po = 4.8. 

Fig. 5. The shape of an initially perturbed regular soliton in the case (T3 = at i = 5, which 
illustrates the stabilization of the soliton through the shedding of small-amplitude radiated waves. The 
plot displayed is the field KeUi{x). The parameters are k = 0.01, a = 1.0,^0 = 0.2, and ui = 0.9. 

Fig. 6. Internal vibrations of an initially-perturbed regular soliton, which was taken close to the 
parameter boundary of the cuspon region. The plot shows the squared amplitude a = \Ui{x — Q)\^ oi the 
Ui{x) field versus time. The parameters are k = 0.01, a — 1.9, t^o = 1-5, and lo = 0.5, with po ~ 0.095 
(see Eq. (H). 

Fig. 7. An example of a stable peakon. The plot shows the field \U\\ versus x and t. The parameters 
are k, = 1.0, a = 1.95, o^o = 1.5, and uj = 0.5, with po = 0.04875. 

Fig. 8. An example of a stable cuspon. The plot shows the field \Ui\^ versus x and t. The parameters 
are k — 1.0, a ~ 1.0, cjo = 1-5, and lu — 0.5, with po = 0.5. 

Fig. 9. The spatial profile of the stable cuspon at t ~ 20. The parameters are the same as in Fig. 8. 
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Fig. 10. The spatial profile is shown for an unstable cuspon in terms of Im?7i at t = 10~^. The inset 
depicts the time-evolution of the maximum value of \Uif. The transition to collapse is clearly seen as 
an explosive temporal behavior of the amplitude. The parameters are k = 0.01, a = 1.1, wo = 0.1, and 
(J = -0.3, with po = 2.618. 
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